A bosonic RVB description of doped antiferromagnet 
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We propose a theory for doped antiferromagnet based on a bosonic resonating-valence-bond (RVB) 
state with incorporating the phase string effect. Both antiferromagnetic (AF) and superconducting 
phase transitions occur naturally within such a bosonic RVB phase. Two distinct metallic regions - 
underdoping and optimum-doping - are also found to be a logic consequence, whose unique features 
explain the recent neutron-scattering measurements in cuprates. 



71.27. +a, 74.20.Mn, 74.72.-h 

It is well-known that bosonic RVB state has a un- 
rivaled high precision in describing the t — J model at 
no-hole case. A mean-field version || of such a bosonic 
spin description also provides a very convenient mathe- 
matical framework. Unfortunately, one soon falls into a 
stalemate once trying to get into metallic phase — always 
encountering a spiral instability jJJ^] at finite doping. 

Recently-revealed phase string effect || has shed new 
light on this dilemma. When a doped hole moves on 
the spin background, the RVB spin pairs are broken up 
along the way and then re-paired, thanks to the superex- 
change coupling. However, such a repairing can never be 
complete for a quantum spin-1/2 media. A phase string 
is found |(| to be left on the hole-path which is nonre- 
pairable in low energy states. It indicates that a nontriv- 
ial Berry's phase will be acquired by a doped hole which 
slowly completes a closed-loop motion. Obviously, such 
a topological effect will be lost if the phase string effect is 
averaged out locally at each step of hopping, which sug- 
gests the spiral phase be an artifact of the inappropriate 
mean- field treatment. 

The phase string effect hidden in the bosonic RVB de- 
scription implies the conventional slave-fermion formal- 
ism H of the t — J model is not a proper starting-point. 
A new exact reformulation || of the model has been ob- 
tained after using a canonical transformation to turn the 
local phase string effect into an explicit nonlocal topo- 
logical effect. Such a new representation is believed to 
be more suitable for a mean-field treatment in both ID 
and 2D cases ||. In this paper, we construct a gener- 
alized bosonic RVB mean-field theory based on this new 
formalism. Main results are summarized in the phase 
diagram shown in Fig. 1, where A s ^ denotes the 
bosonic RVB phase which practically covers the whole ex- 
perimentally interested temperature and doping regime. 
Two real phase transitions occur at low temperature on 
this bosonic RVB background, which are the insulating 
antiferromagnetic long range order (AFLRO) and super- 
conducting condensation (SC). The former is due to a 
Bose condensation (BC) of bosonic spinons and the lat- 
ter is due to a BC of bosonic holons. A striking feature is 
that the AFLRO cannot exist in the metallic phase but 
the spinon BC does persist into an underdoped metallic 



region, leading to a microscopic charge inhomogeneity. 
All of these, including the normal state, are controlled by 
the single bosonic RVB order parameter A s , indicating 
the bosonic RVB pairing to be a crucial driving force re- 
sponsible for a whole spectrum of anomalies in the doped 
antiferromagnet. Some peculiar predications of the the- 
ory will be compared with the experiments in cuprate 
superconductors. 

The t — J model, H t -j = H t + Hj, may be generally 
expressed as 
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In the Schwinger-bonson, slave-fermion representation 
§, one has H i3 = fjfj, B j% = £ CT ufct^, and = 
bi^bj-a- Here is a fermionic "holon" operator 
and bi a is known as the Schwinger-boson operator. The 
bosonic RVB order parameter first introduced by Arovas 
and Auerbach [HI is defined as follows 



(3) 



The sign a appearing in Bji is the source of phase string 
effect mentioned at the beginning. In order to gain a fi- 
nite hopping integral, i.e., (Bji) ^ 0, at the mean-field 
level, we see that up and down spins have to contribute 
differently to avoid cancellation, which always leads to a 
spiral twist. According to previous discussion ||, how- 
ever, this mean-field procedure in doped case may be 
fundamentally flawed as the nontrivial topological ef- 
fect of phase string is totally lost. To avoid this dif- 
ficulty, by a canonical transformation ||, one may re- 
formulate the model such that Bji — J2 a ^"^bj^bia, 
where the singular sign a is replaced by a link variable 
e a*. Correspondingly, Ay and are redefined by 

Hj = Ea e ~^ A ' :3 h*b]-° and H ij = e^hjhj, respec- 
tively. Here the nonlocal gauge field is defined by a 
gauge-invariant condition J2c ^ij = 71 '^c f° r an oriented 
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closed-path C with Nq being the total hole number en- 
closed by C. Apparently vanishes at zero-hole limit. 

And A{ L satisfies J2c A ij =^Y.a° E; e c n t ~ *c with 
n b la = b^bicr and <I>c referring to a uniform flux with 
a strength of 7r-per-plaquette enclosed by C. Here new 
spinon and holon operators, bi a and hi, are both bosonic 
as a peculiar result of the phase string effect. 

Now the phase string effect is precisely tracked through 
the link variables, e icrA ^ and e lA ^ , in B, A s , and H. We 
will still consider the mean-field solution characterized 
by the RVB order parameter defined in Eq. (||). Similar 
to the no- hole case || , the spin part will be diagonalized 
by a Bogolubov transformation 

ha = ^ ( u m1ma - ^mlL-a) e laXm w mty (i) . (4) 
m 

Here "f ma - is an annihilation operator of spinon excita- 
tions and the "single-particle" eigenfunction w mt7 (i) is 
determined by 

j=nn(i) 

with J s — JA S /2. We explicitly introduce a phase factor 
e l,JXm in Eq. (Q) to show a phase uncertainty in w ma 
which cannot be determined by Eq.(^) as it is a linear 
equation. It means that the relation between b and 7 is 
not unique. Without changing A s , such a "phase" free- 
dom can only be fixed by optimizing the hopping inte- 
gral (Bji). Suppose Xm be a general function of the hole 
configuration. Then a maximum hopping (Bji) may be 
achieved by a simple phase shift e i(7Xm — ► — sgn(£ m )e I<7Xm 
each time a hole changes a site 0. The coefficient 
u m and v m in Eq.(@) are given by (\ m /E m + l) 1 / 2 /^ 
and sgn(£ m )(\ m /E m — l) 1 / 2 /\/2, respectively. Here the 
spinon spectrum is E m — yj A 2 ^ — £ ^ , in which the hop- 
ping term only contributes to a shift to the Lagrangian 
multiplier A as X rn = X — Jh/ Js\S,m\- Here the renormal- 
ized coupling constant Jh = (H)t will be always chosen 
as Jh — 5 J below (5 is the doping concentration). A is 
determined by the condition Ei(^L^o-) = N(l — 6), or 

2 -^^L^ coth ^+-BC (6) 

where (3 = 1/ksT and represents the number of 

spinons per site staying at E m = state if a Bose con- 
densation of spinons occurs. In Fig. 1, the region of 
a nonzero A s is shown which is self-consistently deter- 
mined after Ay in Eq. (||) is approximately replaced by a 
mean-field Afy defined by = tt(A^). The effect 

of fluctuations in A;- will be discussed below. 

In the bosonic RVB description, A s does not directly 
correspond to an energy gap, in contrast to the fermionic 
RVB state B. In fact, the spinon spectrum E m is known 



H to be gapless at zero doping and zero temperature 
which ensures a BC of spinons. In the new formulation, 
the transverse spin operator can be written as 0] 

s+ = %M-i)v*< (7) 

in which 

^ = ^Imln (Zi-zi)nl (8) 

describes vortices (with vorticity = 1) centered on holes 
(ri'i = h\hi). In the absence of holes, a Bose conden- 
sation of spinons will always give rise to (Sf) cx (— 1)\ 
i.e., an AFLRO. But in the presence of mobile holons 
— in a metallic region — free vortices introduced by $f 
will generally make (S^) = even though spinons may 
be still Bose-condensed, resembling a disordered phase 
in the Kosterlitz-Thouless transition. Only in an insu- 
lating phase where holes are localized, the AFLRO may 
be still sustained as the vortex effect of $f in Eq. (0) 
can be "screened" through the compensation of a phase 
with opposite vorticities generated from spinons (After 
all, the phase string effect is no longer effective if holes 
are localized). 

So the AFLRO should be absent in the metallic phase 
even though the spinon BC may persist in. One may 
then wonder what is the nature of such a metallic phase. 
In the following, we argue that a spinon BC phase in 
a metallic regime must be generally charge inhomogc- 
neous. Recall that in the BC phase, A takes a value to 
make E m gapless such that n h BC ^ can balance the dif- 
ference between the left and right sides of Eq.(||). Note 
that the E m = state corresponds to the maximum of 
|£ m |, an d thus it is related to those states at the band 
edge of £ m which is generally sensitive to the fluctua- 
tion of A\y As A\j is basically controlled by the holon 
density, the fluctuations of the charge will then leads to 
a "Lifshitz" tail in £ m and play an essential role in de- 
termining the E m — state, which becomes a macro- 
scopic state after the spinon Bose condensation. Such 
a state is generally associated with inhomogeneous hole 
configurations with the condensed spinons forming order 
in hole-deficient region. Obviously, the detailed nature 
of the E m — state will be sensitive to many factors, 
like the dynamics of holons which is beyond the present 
approximation. We will simply treat the fluctuation part 

SA^ = A4J — A? 4 in terms of a random flux descrip- 

l j l j l j 

tion below. Those quantities like transition temperature 
Tbc wm n ot be our primary concern here. Generally 
speaking, with the increase of doping, the reduction of 
the left-hand-side of Eq. ^) will eventually make the BC 
contribution go away. In Fig.l the shaded curve sketches 
such a BC region which basically defines an underdoped 
metallic phase. (The dotted curve in Fig. 1 marks the 
insulating AFLRO phase in the dilute hole regime.) 
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On the other hand, the superconducting order param- 



eter Aff 
following form 



can be expressed M in the 



Af = A^- (fct e 4W-*Jl) (fct e 4fS-*5l) (-1)*, (9) 
where 

$ 4 s =^Imln(z 4 -zOE a < (10) 

describes vortices (anti-vortices) centered on up (down) 
spinons, and <ff( — ► 4h ± 27r after i circles once along 
the loop of a plaquette. In the present bosonic RVB de- 
scription, spinons are always paired (A s ^ 0). In order 
to have a superconducting condensation, then, bosonic 
holons have to undergo a Bose condensation. The vor- 
tices described by $| are all paired up in the ground 
state whose effect is minimal there. But at finite tem- 
perature, free vortices appear in $| as spinons are ther- 
mally excited from the paired state. In order to achieve 
the phase coherence in Eq. ([)]) , the condensed holons have 
to "screen" those free vortices by forming supercurrents. 
A phase transition to normal state eventually happens 
when such "screening" fails which can be estimated as 
the free spinon number exceeds the holon number. The 
BC of holons will be interrupted simultaneously. So the 
transition temperature T c may be determined by 
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K,8, 



(11) 



T=T C 



where k ~ 1 and the left-hand-side represent the number 
of excited spinons determined from Y];„ b ia bj^, T r calcu- 



T c . In turn, the spinon spectrum is qualitatively changed 
which prevents spinons from a Bose condensing into an 
inhomogeneous phase at low-temperature (see below). T c 
shown in Fig. 1 is estimated under such a condition and 
it optimizes T c as compared to the case with stronger 
fluctuations in A^. Beyond T c > Tbc, a crossover due 
to statistics transmutation may quickly set in as holons 
tend to be always Bose-condensed even at high temper- 
ature such that spinons have to be turned into fermions, 
which leads to the breakdown of the bosonic RVB state 
and is beyond the scope of the present paper. 

How can two regions of metallic phase be distinguished 
by experiment? The underdoped metallic region with 
a spinon BC must be a charge inhomogeneity phase 
and also bears some resemblance to "pseudo-gap" phe- 
nomenon. For example, both uniform spin susceptibility 
and transport resistivity will exhibit "pseudo-gap" be- 
havior in this region. But we would like to focus on 
a more direct experimental signature here. The under- 
doped phase may be best characterized by a double-peak 
structure in low-energy region of local dynamic spin sus- 
ceptibility Xl( w ) shown in Fig. 2. Here 

= \ E K ™™ ' t 1 /2(l + n(E m ) + n(E' m )) 

mm' 

-n(E m ,))(u 2 n u 2 n , + vl^,)5{u + E m - E m >)] , (13) 

where u > with K mm > = J2ia \w m a{i)\ 2 \w m >cr{i)\ 2 and 
n(E m ) as the Bose function. The lowest peak in Fig. 2 
originates from excitations from the condensed spinons 
which can be explicitly sorted out from Eq.(13) as 



curve which is obtained under A- 



A% in Eq.©. This 



replacement should be considered to be an "optimum" 
case as the fluctuations of A^ generally reduces T c . We 
will further discuss the optimum condition later. Finally, 
the symmetry of the SC order parameter may be deter- 
mined as 



A? i+ *>/<A? 



i, 



(12) 



where ^ D denotes a summation of A</>9 



5k 



four links of a plaquette and the result indicates a d-wave 
symmetry for the nearest-neighboring SC pairing. 

Therefore, generally there are two temperature scales: 
Tbc and T c , in a metallic phase. At low doping where 
Tbc > T c , the charge inhomogeneity phase happens be- 
low Tbc an d further below T c holons are also expected 
to be condensed into non-uniform regions in favor of the 
spin correlation energy. On the contrary, once T c > Tbc, 
holons will experience Bose-condensation first and be 
uniformly distributed in real space since there is no pre- 
formed local spin ordering. Consequently, Afj may be re- 
placed by Afj with a substantial reduction of 8A\, below 



lated based on Eq.(ll) is plotted in Fig. 1 as the dashed x" ( w ) = ( — n sc ) Z_] Kom——5{u ~ E m ), (14) 

.... . . . . : . \ J. / ^ ' Err 



where m — refers to the BC state. This peak disap- 
pears above TbCi while the second peak is contributed by 
usual spinon pairs excited from the vacuum. Here SA^ is 
treated as a random flux, as noted before, with a strength 
of the flux per plaquette 8<fi — 0.3tt8 at 8 = 1/7 w 0.143. 
Corresponding Tbc is found to be ~ 0.21J. In contrast, 
if the holon BC happens first such that A^- w A^ , then 
the spinon BC is found absent below T c and only a sin- 
gle peak is left as shown in the insert of Fig. 2. Note 
that the sharpness of the peak is due to the Landau- 
level effect caused by A^ in E m (for simplicity we choose 
8(f) = here, i.e., SA^ is totally neglected below T c ). 
Other choices of parameters all yield similar two types 
of structure depending only on whether there is a spinon 
BC or not. 

Neutron-scattering measurement has indeed revealed 
a double-peak structure in Y Ba2Cu 3 Oe ^ 1 compound re- 
cently Q, where the lower peak is located near 30meV 
and the second one is, around 60meV 7 about twice bigger 
in energy as predicted by the theory (Fig. 2). Thus, this 
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underdoped material can be understood in the present 
theory as in the spinon BC phase. On the other hand, 
a "resonance-like" sharp peak at AlmeV has been well- 
known for Y Ba-iCuzO-i below T c pcfl , which is consis- 
tent with the case shown in the insert of Fig. 2 if 
J ~ lOOmeV '. Namely, it corresponds to the uniform 
phase without the spinon BC, defined as an optimum 
metallic region in the theory with an optimized T c . A 
single peak located at Qo =(tt, 7r) is also identified in the 
momentum space at the "resonance" energy in both the 
experiment and theory, suggesting an AF nature of spin 
fluctuations. But an incommensurate momentum struc- 
ture is generally present at energy near the first peak 
shown in Fig. 2 in the spinon BC phase with charge 
inhomogeneity, which is sensitive to the nature of the 
BC state. Finally, it is noted that in the bosonic RVB 
description, the double-layer coupling should not qual- 
itatively change the above energy structure in the odd 
symmetry channel. 



Fig. 1 The phase diagram of doped-antiferromagnet 
based the bosonic RVB description. The dotted and 
shaded curves sketch an insulating AFLRO phase and an 
inhomogeneous metallic region, respectively, described 
by a spinon Bose condensation (BC). SC indicates the 
superconducting condensation region determined under 
an optimal condition (see the text). The temperature T 
is in units of J. 

Fig. 2 Local dynamic spin susceptibility x'Li^) vs - w 
(in units of J) at 8 = 0.143. Solid curve: T = 0; (o): 
T = 0.1; (x): T = 0.2; (*): T = 0.3. Here T BC = 0.21 
with 5cj> = 0.3</>. The insert: x^(w) vs. uj under the 
optimal condition: 8<j> = at T = 0. 
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